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Domains in certain Langmuir monolayers support a texture that is the two-dimensional version of 
the feature known as a boojum. Such a texture has a quantifiable effect on the shape of the domain 
with which it is associated. The most noticeable consequence is a cusp-like feature on the domain 
boundary. We report the results of an experimental and theoretical investigation of the shape of a 
domain in a Langmuir monolayer. A further aspect of the investigation is the study of the shape 
of a "bubble" of gas-like phase in such a monolayer. This structure supports a texture having the 
form of an inverse boojum. The distortion of a bubble resulting from this texture is also studied. 
The correspondence between theory and experiment, while not perfect, indicates that a qualitative 
understanding of the relationship between textures and domain shapes has been achieved. 
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I. INTRODUCTION 

The direct observation of monolayers at the air/water 
interface by the methods of polarized fluorescence mi- 
croscopy (PFM) |l| and Brewster angle microscopy 
(BAM) Q reveals that the films possess complex textures 
similar to those observed in liquid crystals. The textures 
are generally found in "tilted" phases, i.e. in phases in 
which the long axes of the molecules in the film are are 
not perpendicular to the water surface but are uniformly 
tilted with respect to the normal. The textures are the 
result of the spontaneous organization of the molecular 
tilt azimuth on macroscopic length scales. They can be 
understood, at least qualitatively, in terms of continuum 
elastic theories of smectic liquid crystals ||. 

Many striking textures have been found in domains 
of condensed tilted phases, such as the Li phase, that 
are surrounded by an isotropic phase, either liquid [often 
designated LE, (liquid-expanded) or L{\ or gas (G) 
Among those textures are boojums, in which the tilt az- 
imuth varies continuously and appears to radiate in some 
cases from a defect located at the edge of the domain and 
in others from a "virtual" defect in the isotropic phase 
[^j. The shapes of domains that contain the boojum tex- 
ture are not round; they have a cusp-like deformation 
pointing toward the defect. 

In essence, the equilibrium shape of a domain contain- 
ing a boojum is like the equilibrium shape of a crys- 
tal, which can be calculated with the use of the Wulff 
construction Q. In that procedure it is assumed that 
the surface energy of the crystal is determined solely by 
the anisotropy of the bulk energy. For a monolayer do- 
main containing a nontrivial texture this is not the case; 
the boundary — more properly, line — energy is compara- 
ble to the energy associated with the alignment of the 



tilt azimuth, so shape and texture must be calculated 
self-consistently. 

Such a calculation has been performed by Rudnick and 
Bruinsma , who showed that there is a cusp on the do- 
main boundary, quantified by an "excluded" angle "J? that 
varies with R, the radius of the domain. Their results 
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FIG. 1. (a) The excluded angle, St, plotted as a function of 
the 1/R. (b) The definition of the excluded angle 
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are displayed concisely in Fig. [l], which contains a pic- 
torial definition of '5 and a plot of ^ versus 1/R. The 
excluded angle is equal to zero when the domain is cir- 
cular and increases as the cusp sharpens. Rudnick and 
Bruinsma predict that '5 vanishes in the limit R — oo and 
that the angle is initially linear in R -1 . The excluded an- 
gle goes through a maximum and then decreases to zero 
for sufficiently small domains. Schwartz, et al. (§) made 
PFM measurements of ^ for domains of pentadecanoic 
acid ranging in size from 12 to 120 //m. They found that 
the excluded angle increases with 1/R, but roughly as 
a power law with an exponent of 0.3 rather than unity. 
There was no evidence of a maximum. 

Riviere and Meunier || carried out BAM studies of 
condensed domains of hexadecanoic acid surrounded by 
a gaseous phase and measured the distance, d, of the vir- 
tual defect from the domain edge as a function of R. By 
assuming that the deformation of the domain is small, 
they were able to derive a relation between the quan- 
tity p = d + R and the domain radius from continuum 
elastic theory. They compared the calculated and mea- 
sured plots of p versus R and were able to obtain values 
for the ratio of the bend/splay elastic constant to the 
anisotropic part of the line tension. They also made a 
qualitative comparison between the measured elongation 
of the boojum-containing domains and the elongation 
calculated from a generalized Gibbs-Thompson equation 
that was obtained by Galatola and Fournier jl0| , this un- 
der the assumption that the texture is that of an undis- 
torted boojum. 

Experimental tests of the theory have been limited by 
the small range of domain sizes that could be examined 
by optical microscopy. In this paper we report new mea- 
surements of the excluded angle for much smaller do- 
mains. This has been accomplished by transferring the 
monolayers to a solid support by the Langmuir-Blodgett 
technique and obtaining images with the use of scanning 
force microscopy. We also describe measurements on "in- 
verse" boojums, which occur when there is a "bubble" 
of isotropic phase imbedded in a tilted phase. A new 
theoretical analysis of the boojum problem, which is an 
extension and reformulation of the treatment of Rudnick 
and Bruinsma, is then described. Comparisons are made 
between the experiments and this theory. 



II. EXPERIMENT 

Monolayers of pentadecanoic acid (Nu-Chek Prep, > 
99%) were deposited from chloroform (Fisher spectran- 
alyzed) solutions onto water (Millipore Milli-Q, pH 5.5) 
in a NIMA Type 611 trough. They were transferred onto 
mica substrates that had been freshly cleaved with ad- 
hesive tape and immediately inserted into the water sub- 
phase. The transfers were performed by withdrawing the 
mica from the subphase at constant speeds ranging from 
0.5 to 2 mm/min. Boojums of the L 2 phase surrounded 



by the LE phase were obtained in the coexistence region 
at temperatures between 15 and 25°C; boojums of the 
G phase surrounded by the L 2 phase were obtained in 
transfers at 7°C. 

The transferred monolayers were imaged with a scan- 
ning force microscope (SFM) (Park Scientific Instru- 
ments) using a 100-/im scanner at room temperature in 
ambient atmosphere. The instrument was calibrated on 
the micrometer scale by imaging grids with known spac- 
ings and on the nanometer scale with mica standards. A 
microfabrication triangular SiaN4 cantilever with a nor- 
mal spring constant of 0.05 N/m was used for the mea- 
surements. The vertical bending and lateral torsion of 
the cantilever were monitored by reflecting a laser bean 
from the end of the cantilever into a four-segment pho- 
todetector so that the topographic and frictional force 
images of the samples could be obtained simultaneously 
and independently of each other. All images were ob- 
tained in the constant-force mode. The loading forces 
ranged from 1 to 5 /zN. 

Brewster angle microscope images of Langmuir mono- 
layers were obtained with an instrument of our own con- 
struction. Details of the instrument and the experimental 
procedures can be found elsewhere |yj . 



III. EXPERIMENTAL RESULTS 



The key to the success of the experiment is the re- 
quirement that the boojums survive the transfer process 
and that their shape on the support be the same as that 
on the water surface. It has often been observed that 
the continuous LE phase breaks up into small islands 
when it is transferred to a solid support. The behavior 
has been found in monolayers of stearic acid deposited 
on thin polyethyleneimine films [fl2f , dimyristolyphospha- 
tidic acid deposited on glas s fili and octadecyltrichlorosi- 
lane on acid-treated mica |14| . However, as can be seen 
in Fig. ||(a), despite the breakup of the LE phase, the 
boojum shapes of the L 2 phase survive. LB transfer is 
also known to cause flow in monolayers |l5|,[l6| , and there 
is evidence of this in the transfers made at 2 mm/min. 
As is evident in Fig. ||(b), the domains are elongated in 
the dipping direction. However, at 0.5 mm/min no dis- 
tortion is evident and the boojums are not aligned in the 
dipping dirccion. 

A plot of the excluded angle against the reciprocal of 
the domain diameter is shown in Fig. J3J. The same data 
are displayed in a log-log plot in Fig. [J In the range 10 
- 20 /im the data from the SFM studies agree with those 
obtained earlier by PFM ||. We note that the slope of 
the log-log plot at large R is less than unity, and that 
there is a non-zero intercept in the '5 vs 1/R plot. The 
fact that the intercept is not zero had been overlooked in 
the earlier work, in which the data were presented on a 
log-log plot. By extending the measurements to smaller 
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FIG. 2. Atomic force microscope images of domains asso- 
ciated with boojums. Monolayers were transferred onto mica 
from the L2 — LE coexistence region at a speed of 0.5mm/min 
(a), (c), and at 2 mm/min (b). 

domains we have observed that there is a maximum and 
that the small domains become circular. There is no evi- 
dence of a cusp in domains smaller than 2 ^m in diameter; 
see Fig. ||(c). 

The good agreement between the PFM and SFM stud- 
ies for the larger domains gives us confidence that the 
shapes have not been seriously changed during the trans- 
fer process or greatly affected by the change in substrate. 
We cannot be certain that this is true of the smaller do- 
mains, but it does seem likely that the effect of hydrody- 
namic flows and variations in the substrate will be most 
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FIG. 3. The excluded angle of the cusp-like feature on a 
domain supporting a boojum texture, plotted against the re- 
ciprocal of the domain radius. Filled triangles represent data 
for domains on water surface and empty squares represent 
data for domains tranferred onto solid support. 
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FIG. 4. A log-log plot of the dependence on radius of the 
excluded angle of the cusp-like feature on a domain support- 
ing a boojum texture. Filled triangles represent data for do- 
mains on water surface and empty squares represent data for 
domains tranferred onto solid support. 

serious in the larger domains. 

In the PFM studies Q it was observed that inclu- 
sions of the LE phase in the L2 domains also had cusps. 
When we examined films at temperatures below the LE- 
L 2 -G triple point at 13° C [0, we found similar cusp- 
shaped regions in the isotropic G phase surrounded by 
the anisotropic L2 phase. These domains also survive 
transfer to a solid support. Fig. [| is a frictional force 
image of such an "inverse" boojum. The frictional force 
within the "inverse" boojum is high because the tip im- 
ages the mica substrate and not the dilute G phase. This 
is also confirmed by the topographic image, which was 
obtained simultaneously. The height profile shows that 
the inverse boojum extends to the substrate surface. 
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IV. THE INVERSE BOOJUM 



S |lttl 

FIG. 5. Frictional force image of an inverse boojum consist- 
ing of a bubble of the isotropic G phase surrounded by the 
anisotropic L2 phase. The monolayer was transferred onto 
mica from the L2 — G coexistence region at a speed of 0.5 
mm/min. 
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FIG. 6. The excluded angle of the cusps in bubbles of 
isotropic G phase surrounded by anisotropic L2 phase, plotted 
versus the reciprocal of the bubble radius, filled triangles rep- 
resent data for bubbles on water surface and empty squares 
represent data for bubbles tranferred onto solid support. Note 
the similarity with Fig. 3. 

Measurements of the excluded angles in these inverse 
boojums as a function of their radii are plotted in Fig. 
^. The behavior is similar to that found for the ordinary 
boojums; the excluded angle has a nonzero intercept, in- 
creases roughly as R~ 3 , goes through a maximum and 
falls to zero. Included in the figure are three points ob- 
tained from BAM studies of inverse boojums in Langmuir 
monolayers of pentadecanoic acid. There is a reasonable 
agreement with the SFM experiments for large domains. 
However, cusp angles are difficult to measure in small 
domains by BAM. They can be more easily and precisely 
measured by fluorescence microscopy, but the precipita- 
tion of the probe in the L2-G two-phase region made this 
impossible. 



The boojum texture is an outcome of the competition 
between bulk and boundary energies in a bounded do- 
main containing a condensed phase. When the boundary 
energy has the particular anisotropic form displayed be- 
low, and the domain is a perfect circle, this texture is an 
exact solution to the energy minimization equations [M. 
At the same time as the XY-like order parameter adjusts 
in response to the imperatives of energy minimization, 
the boundary also distorts so as to accomodate the now 
inhomogeneous environment presented by the texture. It 
has been f oun d that when the energy is as given by Eqs. 
(4.1) and (4.4), a circular domain satisfies the extremum 
equations. 

This means that the simple version of the energy of 
a domain in a Langmuir monolayer leads to a liquid- 
condensed-phase domain supporting a nontrivial texture 
whose boundary is perfectly circular. To explain the do- 
main shapes observed in the experiments reported here 
and previously ^,^| it is necessary to introduce elabora- 
tions on the system's energy. The addition of a higher- 
order harmonic term to the line tension does give rise to 
a distortion of the boundary and a feature in the form 
of a cusp appears. The theory is in qualitative accord 
with the experiments except in the limiting large-domain 
behavior, where it predicts that the cusp angle vanishes 
rather than going to a limiting non-zero value as observed 
for both domains and bubbles. 

The shape of bubbles of an isotropic phase inside an 
ordered phase was not considered in the theory. When 
there is a two-dimensional bubble of gas phase inside an 
extended region of the liquid-condensed phase a texture 
appears in the liquid-condensed phase that is closely re- 
lated to the virtual boojum observed in two-dimensional 
liquid-condensed droplets. This texture minimizes the 
total energy H [Q(x, y)], given by 
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e(x,y))ds, (4.1) 
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where Q(x,y) is the angle between the vector order pa- 
rameter c and the x-axis, while 9 is the corresponding an- 
gle for the unit normal, n. The coefficient n is the Frank, 
or "spin stiffness" constant. The energy above follows 
from the assumption of equal bend and splay moduli, 
which is known to represent, at best, an approximation 
to physical reality. The function a (6 — Q(x,y)) is the 
surface energy, which depends on the relative orienta- 
tions of c and the unit normal to the boundary, h. The 
extremum equations that arise are 



V 2 e(x,y) =0 

on 



(4.2) 
(4.3) 



Eq.(4.2) applies outside the bubble and Eq.(4.3) at the 
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FIG. 7. The texture associated with an inverse boojum and 
t he p osition of the virtual singularity, given by Eqs. (4.5) and 
(^q). The distortion of the bubble from a circle, obtained 





from Eq. (5.24), is also indicated on the Figure. The param- 
eters used are k/ol = 0.4/xm, a/erg = 0.16 and R = 2.4/xm. 



interface. If cr(</>) has the form 

<j(4>) = uq + a cos 4>, 



(4.4) 



and the bubble is a perfect circle, then a solutio n fo r 
Q(x,y) t hat satisfies the extremum equations Eq. (4.2) 
and Eq. (JO) is 



Q(x,y) 



In 1 



x — ty 



In 1 



x + ly 



(4.5) 



Given the similarity between the mathematical structure 
of the right-hand side of Eq. (4.5) and the expression for 
Q(x, y) in the case of the virtual boojum (see ref. 0) this 
solution can be characterized as an inverse boojum. It is 
a texture with a singularity that lies a distance Rb from 
the origin, where 



Rb 
Ro 



a 

Ro 



clRo/k 



yf 1 + {aRo/nf + 1 



(4.6) 



The above relation serves to fix the value of the quan- 
tity a. Note that Rb is always less than the radius of 
the bubble. The singularity lies within the bubble. The 
texture associated with a virtual boojum is displayed in 
Fig. 

When a virtual boojum exists in a condensed-phase do- 
main, the orientation of the director leads to an optical 
anisotropy that should be observable by Brewster-angle 
microscopy pf)). A "sunburst" of darker and lighter 




FIG. 8. The curves along which the order parameter points 
in a fixed direction when there is an inverse boojum texture 
in the immediate vicinity of a bubble. The parameters used 
are n/a = OAfim, a/ao = 0.16 and Ro = 2.4/im. 

lines radiates from the point of singularity. Regions of 
constant intensity on the pattern correspond to portions 
of the texture in which the order parameter points in a 
fixed direction. In the case of the inverse virtual boojum, 
there should be an analogous signature of the t extu re in 
the surrounding condensed phase. Using Eq. (fT^) one 
finds straightforwardly that the orientation of the order 
parameter is constant along the perimeter of a circle. 
This means that points of equal orientation fall on circu- 
lar countours as shown in Fig. |^. The contours appear 
to pass through two points inside the bubble. One of the 
points is the location of the singularity, and the other 
is the center of the bubble. The optical anisotropy that 
will be observed by BAM depends on the molecular pa- 
rameters and the orientation of the boojum with respect 
to the direction of the incident light. Qualitative simu- 
lations of the images are shown in Fig. ^[ They have 
been obtained by using the relation between tilt azimuth 
and reflectivity calculated by Tsao, et. al. Jig] . These 
images can be compared with the BAM images in Fig. 
|To| , which show similar variations in the reflectivity. It is 
interesting to note that the boojums appear to align. A 
similar alignment of isotropic droplets in a nematically 
ordered host phase has been observed Jl9|]. This align- 
ment is attributed to a long-range interaction produced 
by distortions in the director fields. 



V. BUBBLE SHAPE: THEORY 



In equilibrium a bubble takes on the shape that mini- 
mizes the system's free energy. If the interior is uniform 
and the line tension of an element of boundary depends 
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FIG. 9. Simulated BAM images of inverse boojums for dif- 
ferent orientations with respect to the incident light. The 
orientation of the beam with respect to the symmetry axis of 
the boojum is specified in the lower corner of each figure. 

on the orientation of the boundary element, the Wulff 
construction can be utilized to yield the shape. In the 
case of inhomogenous interior structure, a generalized 
version of the Wulff construction provides an approach 
to the solution J?j] . This method applies to the bubble in 
much the same way as it does to a domain. One discovers 
immediately that a circular bubb le does not minimize the 
system energy of the form ( fl.l[ ) with <r((f>) given by Eq. 
( [4.4[ ), in contrast to the result obtained for a domain 0, 
for which a circle represents the minimum energy domain 
shape. However, the method cannot be applied directly 
because the curvature of the boundary changes sign. As 
a consequence, the radius of curvature of the boundary 
may become infinite at points along it, and this proves 
awkward in the context of the Wulff construction. How- 
ever, a reparameterization allows one to proceed with 
the analysis. We find that there is no feature that can, 
strictly speaking, be described as a cusp, in that the slope 
of the bounding curve has no discontinuity. There is, on 
the other hand, a feature on the boundary whose shape 
can be quantified in terms of an excluded angle. 



We have utilized an alternative approach to the deter- 
mination of the bubble's shape, to provide a check on 
our Wulff construction analysis. We adopt polar coordi- 
nates. This formulation confirms the results of the anal- 
ysis based on the Wulff construction. We obtain explicit 
results for the bubble's shape. Cusp angles are deter- 
mined both analytically and by means of measurements 
on figures generated with the use of explicit analytical re- 
sults, and these angles provide the basis for comparison 
with experiment. 



FIG. 10. BAM images of inverse boojums in Langmuir 
monolayers of pentadecanoic acid at 5°C. 

The remainder of this section consists of a review of 
the generalized Wulff construction, the development of 
the alternative approach and a summary of results. De- 
tails are deferred to a forthcoming publication p0| . 



A. Generalized Wulff construction 

Following the approach of Rudnick and Bruinsma, we 
attempt to apply the generalized Wulff construction Q] 
to determine the shape of the bubble. In our applica- 
tion, this technique is used to minimize the total energy 
Eq. (4.1) with respect to the Burton, Cabrera and Frank 
(BCF) parameterization of bubble shape, i.e. R{9) and 
9 |§. 

We rederive the Wulff construction so as to remove 
the approximation — utilized in earlier work Q| — that the 
radius of curvature, R + R" , is constant. We obtain 
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It can be shown that if <j(4>) — a(—(p), we have for the 
bubble 
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a' (cost/ — — h suit/— ) 
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Adopting (4.4) and (4.5), the equation for the shape 
of the bubble is given by 



+cr (- sin 9— + cos fl- 
ax 



dy' 



9, 



(5.2) 
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2a(a + S) cos 29 - 4a(l + a 2 + 25a) cos 9 + 1 + 4a 2 + a 4 + 5a(l + 6a 2 - a 4 ) ' 



(5.3) 



where 5 — a/ag, R = Rqt and i?o is the ra dius of the 
unperturbed circle. The denominator of Eq. ( |5.3| ) passes 
through zero in regimes of interest, and this complicates 
the analysis. It proves useful to reparameterize the sys- 
tem in Cartesian coordinates. 

The relations between the parameterization using 
Cartesian coordinates and the BCF parameterization are 
as follows: 
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Using Eq. ( |5.4| ) and Eq. ( p.5| ), we have the following 
relations for nearly circular boundaries and in the region 
where 9 is small 
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obtain 
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The integration of Eq. ( |5.§| ) is readily performed easily 
to yield 
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and 
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The next step is to fix the constants of integration 
K y i and K y . Symmetry requires y'(—x) = —y'(x), and 
yields K y > — 0. The bounding curve y(x) will merge into 
a circle of radius 1, centered at the origin. The merg- 
ing point, (xq, y(xo)), is taken to be the point at which 
y'(xo) — —xqI\/\ — x\. This leads to 



Ky = 

C 2 x 2 2(1 -af + da 
2a 2 + 2a 3 



ln[aa% + (1 - a) 2 } 
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2(l-a) 3 al 
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+ -\/l-x 2 , 



a 2 xo 
1 — a 

(5.15) 



The boundary of the bubble Eq. ( |5.14 ) is smooth in 
that there is no discontinuity in the slope of the bound- 
ing curve. For very small Rq, the radius of curvature of 
the boundary close to the virtual boojum, b = l/y"(0), 
is approximately equal to i?o- The bubble appears to 
be circular and the excluded angle, ^, is zero. The ra- 
tio (5 = b/Ro decreases as Rq increases. This causes the 
bubble to sharpen at one end. When j3 exceeds a critical 
value 0c, this sharper end has properties in common with 
be a cusp, and one is able to derive a excluded angle, 
as depicted in Fig. |TT|a. It can be shown that \t is given 
by the following: 



* = 2 cos" 



h ' 
1-6 



(5.16) 
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where h — \y(0) + l|is the deviation of the tip of the 
feature from a circle. Because there never a cusp in the 
strict sense of the term the choice of (3 C , is a matter of in- 
dividual judgement. Based on theoretical bubble shapes 
obtained in Fig. [l^, f3 c is chosen to be 0.45. 

When i?o ^> k/cl, the x 2 term in Eq. 
nates. 



(5.13) domi- 

The boundary around the cusp-like feature is this 



x term instead of the unperturbed circle. The new def- 
inition of the excluded angle is as displayed in Fig. pd| b 
and we find W can be expressed as 




m(l - a) 



(5.17) 



In this limit, a — > 1 — e, where e oc r 1 < 1, we approxi- 
mate C\ as follows 



Ci «2e(e-2tf) + 0(5e 2 ), 



(5.18) 



We then have the large- i?o behaviour of the excluded an- 
gle 



In - 
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FIG. 11. (a) The way in which 9 is defined where the 
boundary is smooth. This method leads to an analytical ex- 
pression for the excluded angle of a cusp-like feature, (b) For 
Ro 2> n/a, the "cusp"-like feature show up differently and 
requires a new definition of 9. 



B. Polar Coordinates 

As an alternative, one can parameterize the domain 
boundary in polar coordinates. In this case, we choose an 
origin, nominally at the center of the undistorted circular 
bubble. The cartesian coordinates x and y are replaced 
by r, the distance from the origin and 4>, the angle with 
respect to the x axis. In terms of r and cj>, the inverse 
boojum texture has the form 



Q(r,<f>) = \ 



Infl-^-lnfl-^-* 



(5.20) 



The boundary condition, Eq. (T3) leads to the following 
result for the parameter a: 



a 



1 



^{n/aRof + 1 + (n/aRo) 



(5.21) 



where Rq is the radius of the circular bubble. The param- 
eter a has a simple interpretation — it is the distance of 
the singularity associated with the virtual boojum from 
t he c enter of the bubble. Note that according to Eq. 
(|5.21 ) the singularity always lies within the bubble. 

We now proceed with a calculation of the distortion 
of the bubble from a circle. Our results are the leading 
order terms in an expansion in the dimensionless, and 
presumably small, combination a/a . If one writes 



r(4) =R e k ^\ 



(5.22) 



one eventually arrives at 
R a k"(<p) = 



1 



1 



-i<p 



(5.23) 



Eq. (5.23) represents the leading term in an expan- 
sion in the dimentionless combination a/do- This ra- 
tio is small, as the line tension is assumed to be nearly 
isotropic. 



The double integration of Eq. (5.23) is relatively 
straightforward. One finds for k((f>): 



k(4>) 




arctan 



Ro 



1 — COS ( 

Hn 



+2 In 1 



a \ a 
TT -2— cos< 

Kq I -rto 



a 
Rq 



(5.24) 
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(a) 




(b) 




FIG. 12. (a) The calculated shapes for the bubble of 
7?o=0.4, 0.8, 1, 1.2, 1.6, 2, 2.4 and 3.2 fj,m. The bubble 
which is drawn in thicker line is the one of the nominal size, 
R c = lfim when /3 C = 0.45 for this case, (b) The calculated 
shapes for the bubbles of 7?o=3.2, 8 and 16 fim. The param- 
eters used are n/a = 0.4/xm and a/ero = 0.16. 

In the regimes in which comparisons can be made, the 
results of the calculation based on the utilization of polar 
coordinates are in complete accord with results that fol- 
low from the Wulff construction. In addition, an analy sis 
based on polar coordinates — in particular Eq.( 5.24| ) — 
provides figures that can be utilized for comparison with 
experim ental data. Bub bles w hose shapes are described 
by Eq. ( |5.22j ) and Eq. ( |5.24|) are displayed in Fig. [if. 
These bubbles comprise the results of theoretical model- 
ing that are compared with the data described in earlier 
sections. 



VI. RESULTS AND CONCLUSIONS 



We now compare the theoretical excluded angles with 
the experi menta l ones. Within the range of experiment 
data, Eq. (5.16) provides the basis for comparison be- 
tween experiments and theory. We find that the best fit is 
achieved when the parameters are chosen as n/a — 0.4/im 
and a/do = 0.16. The comparison between theory and 
experiment is displayed in Fig. [l3[ Because of the fact 
that there is never a cusp in the form of a discontinuity 
in the slope of the bounding curve, we are unable to 



60.0 



40.0 




20.0 



0.4 0.6 0.8 

1/R(um) 

FIG. 13. The experimental measured excluded angle of the 
cusps in bubbles, as in Fig. 6, shown as empty squares and 
filled triangles, the measured excluded angle of the calculated 
bubbles shown as empty diamond and the analytic expression 
of the theoretical excluded angle shown as solid line, plotted 
vesus the reciprocal of the bubble radius. Note the fundamen- 
tal difference between the behavior of the the small-radius be- 
havior of the theoretical curve, which tends monotonically to 
a finite limit at vanishing bubble radius and the measured ex- 
cluded angle of calculated bubbles, which seems to rise rapidly 
from zero as the bubble radius exceeds a "threshold" value. 

definitively model the "onset" of a cusp. However, we do 
find an apparent rapid dropoff to zero of the (nominal) 
excluded angle. Furthermore, the theory predicts that 
sufficiently small bubbles will be nearly circular. It is 
possible to determine a value of what appears to be the 
radius of onset, i? c , with the exercise of a little judge- 
ment. For the parameters of choice, we find R c = Ifim 
for p c = 0.45. 

Our theory appears to reproduce fairly well both the 
qualitative and the quantitative behaviour of the bubble 
shapes as a function of the bubble size. Small bubbles 
are nearly circular and the excluded angle, it, is equal 
to zero. Cusp-like features are observed for bubbles hav- 
ing radii Rq greater than a threshold value R c . The an- 
gle 'J decreases for bubbles of larger sizes. However, in 
the experiments, the excluded angle '5 tends to a non- 
zero constant as Ro — > oc. This limiting behavior is 
not achievable within the framework of the theory pre- 
sented here. It is known that thermal fluctuations will 
cause the line tension coefficients to depend on the diam- 
eter of the bubble or domain p3|] . However, preliminary 
calculations indicate that this renormalization has no ef- 
fect on the asymptotic behavior of the excluded angle 
for large domains of bubbles. A possible source of the 
discrepancy is the difference between the bend and splay 
moduli, which has been neglected in the present work. 

The ratio n/a for hexadecanoic acid has been obtained 
by Riviere and Meunier M. They find 28/xm < n/a < 
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40//m, about two orders of magnitude larger from our 
best estimate, n/a — OAfim. Studies of other textures in 
monolayers such as stripes |p3 indicate that k increases 
with chain length. This could account for some of the 
difference. Another possible reason for this discrepancy 
is the fact that Riviere and Meunier adopted n\ ^ K3 as 
the basis for texture distortions while we have used a, 
the anisotropic line tension, as the mechanism for bub- 
ble shape distortions. Our comparison between theory 
and experiment gives an estimate of a/ao = 0.16 for 
pentadecanoic acid at LC-G boundary. The line ten- 
sion anisotropy, which corresponds to of our a/do, has 
been estimated for the LC-LE interface in D-myristoyl 
alanine from studies of dendritic growth ^5|. The value 
found is 0.005 - 0.025. Given the fact that the systems 
are not identical, and that these are only estimates of the 
anisotropy, the difference is not at all disturbing. 

In summary, we have formulated a theory for bubbles 
in Langmuir monolayers that is an extension and refor- 
mulation of the theory of domain shapes discussed by 
Rudnick and Bruinsma 0. Although the current theory 
predicts a smooth bubble boundary, there exist cusp-like 
features with well-defined and measurable excluded an- 
gles when the radius of the bubble exceeds a nominal 
value R c . Comparing theoretical excluded angles with 
the results of experiments, we find good qualitative and 
quantitative agreement. However, there are also discrep- 
ancies. Experimental observations and related work indi- 
cate that the difference between k± and K3 cannot be ne- 
glected. The calculation of an equilibrium texture when 
K\ 7^ K3 is currently being performed as a prelude to the 
determination of the shape of a bubble or domain. 
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